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Abstract 

T-|- ' Results for higher-order threshold enhancements in high-E'r jet production in hadron- 

CN . hadron collisions are presented. Expressions are given for the next-to-next-to- leading order 

(NNLO) threshold corrections to the single-jet inclusive cross section at next-to-leading 
^ ' logarithmic (NLL) accuracy. The corrections are found to be small for the specific choice 

OO . of EtJI for the factorization and renormalization scales, and the corrected cross section 

shows a substantial reduction of the scale dependence. A comparison to experimental 
[*-. ' results from the Tevatron is presented. 
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1 Introduction 

QCD provides an impressive description of the data for the production of high-i^^ jets in 
hadronic colhsions. The data taken by the CDF [|I| and D0 Collaborations for the single jet 
inclusive cross section fall by nine orders of magnitude in the Ex range extending to 460 GeV, 
a behavior also seen in the theoretical predictions. Nevertheless, there is a hint that the theory 
slightly underestimates the data towards the upper end of the measured Et range. Although 
such an excess of experiment over theory could be a sign of new physics, other more conventional 
explanations are possible. On the one hand, it could be that the parton distributions at the large 
values of x needed for the theoretical calculations are slightly underestimated. In particular, 
a somewhat harder gluon distribution than is usually used might be sufficient to describe the 
data. This approach has been studied by the CTEQ Collaboration [Q which developed a special 
parton set, CTEQ5HJ [^], to describe these data. On the other hand, it may be that there are 
corrections to the hard scattering parton-parton cross sections that have yet to be included in 
conventional next-to-leading-order calculations. It is this latter possibility which is examined 
in this paper. 

At large values of xt = 2Et/vS, the phase space for gluon emission is limited by the rapid 
decrease of the parton distributions. The parton-parton scattering processes are restricted to 
the threshold region where there are large logarithmic corrections. Threshold resummation is 
a formalism for including the effects of these corrections to all orders in perturbation theory. 
The effects due to threshold resummation are expected to grow as xt approaches one. The 
high-i?r jet data extend in xt to about 0.5 and it is of interest to see if some of the observed 
excess of data over theory can be accounted for by the inclusion of threshold corrections. 

In the next section, the formalism for threshold resummation to NLL accuracy for single jet 
production in hadron-hadron collisions is reviewed. Sec. 3 contains a comparison of our results 
with Tevatron data while Sec. 4 contains our conclusions. Detailed expressions for the NNLO 
threshold corrections, derived from the expansion of the resummed cross section, are given in 
a set of appendices for each of the relevant parton-parton subrocesses. 

2 The single-jet inclusive resummed cross section 

In this section we discuss the single-jet inclusive cross section at large transverse momentum, 
paying particular attention to threshold corrections that occur at each order in perturbation 
theory. Consider jet production in pp collisions, 

p + p^J + X, (2.1) 

the invariant single-jet inclusive cross section for which may be written in f actor ized form as 

X Ej^-^^^^^^is, t, M, /ip, asiiJ.1)) ■ (2.2) 

d^pj 

The initial-state collinear singularities have been absorbed into the parton distribution functions 
at a factorization scale ^f, while ^r is the renormalization scale. The parton-parton hard 
scattering cross section is denoted by a. 
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The parton-parton scattering subprocesses contributing to jet production in lowest order 
are 

mi , Ijik -^ QjQk , 

(2.3) 

For the process faiPa) + fhiph) -^ Jipj) + ^, the Mandelstam invariants constructed from the 
parton and jet four- vectors are given by 

S={,Pa+Phf , t={.Pa-Pjf , U={pb-Pjf , (2.4) 

which satisfy S4 = s + t + u = 0at threshold. The variable s^ is the square of the invariant 
mass of the system recoiling against the observed jet. 

The threshold for the above partonic subprocesses occurs at s^ = 0. The values of Xa and Xf, 
corresponding to this point are the minimum values which can give rise to a jet of the specified 
rapidity and transverse momentum. Near threshold the phase space for the emission of real 
gluons is limited and large logarithmic corrections arise from the incomplete cancellation of 
infrared divergences against virtual gluon emission contributions. In general, a includes plus 
distributions with respect to S4 at nth order in Og of the type 



S4 



m < 2n - 1 , (2.5) 



where pt = {tu/sY^"^ is the transverse momentum of the jet. These distributions have been 
resummed to all orders at next-to-leading logarithmic (NLL) accuracy for dijet and singlet-jet 
production in Refs. [^ |, 0. 

The resummation is achieved in moment space through a refactorization [^ |3, H, of the 
cross section into a product of functions ip and J that absorb the collinear singularities in the 
incoming partons and outgoing jets, respectively; a soft gluon function S that encompasses 
noncollinear soft gluon emission; and a hard-scattering function H that describes the short- 
distance hard-scattering. The color structure of the hard scattering depends on the flavor 
content (for example singlet or octet for qq -^ qq) and is described by a set of color tensors 
cj. Note that the functions S and H are actually matrices in color space. At lowest order 
Sli = Tr[c\^ci]. As we shall see, the trace of the product of the matrices H and S at lowest 
order gives the Born cross section. Full details of the resummation formalism for jet production 
are given in Refs. §, |, 0, |TD|, [TTl. 

The moments of the perturbative partonic cross section are then given by 



(Ppj J s cP'pj 



= ^^°/^-i^"|^^-/^-|^\^ J»>(iV)/7g(p^/(Ar^^)) . (2.6) 

where J represents the final-state observed jet, and J*" the partons recoiling against the jet. 
Note that Na = N{—u/s) and Nb = N{—t/s), with A^ the moment variable. The moments 
of the plus distributions, Eq. (|2.5| ), in the cross section produce powers of InA^ as high as 
In " A^. The leading and next-to-leading logarithms of A^ are then resummed to all orders in 
perturbation theory. 



The resummed cross section in moment space follows from the resummation of the A^- 
dependence of each of the functions in Eq. ( |2.6|) , and may be written as f^, 0] 
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The first exponent resums the iV-dependence of the incoming partons [jT2|, \13 
leading order accuracy in InA^, A^^^ag) = C/[(«s/7r) + (_ft'/2)(as/7r)^], and u^^^ 
Here Cf = Cp = {N^ — l)/(2A'"c) for an incoming quark, and Cf = Ca = Nc for an incoming 
gluon, with Nc the number of colors, while K = Ca(67/18 — 7r^/6) — 5nf/9, where n/ is the 
number of quark flavors. Also C^ = p'j/pr- Note that the leading contributions from this 
exponent enhance the cross section [0 . 

The second exponent resums the A^-dependence of the final-state observed jet, J, and of par- 
tons r recoiling against the jet [|,|^, IH, with B[^^ = (a,/7r)(-3Cp/4), BJ. = (a,/7r)(:7p[ln(2z/g)- 
1], B[g^ = («,/7r)(-/5o/4), and B'^^^ = (a,/7r)C^[ln(2z/g) - 1] 0, |T|, |ll]. The Ui terms are 

gauge dependent. They are defined by Vi = {Pi ■ n)'^/\n\'^, where Pi = piJ2/s are the par- 
ticle velocities and n is the axial gauge vector, chosen so that Pi ■ ( = Pi ■ n ioi i = a,b 
0, |l5l and |rap = 1. Also Pq = (IIC^ — 2nf)/3 is the one-loop coefficient of the P function, 
P{as) = fidlng/dfi = — /5oas/(47r) + .... Note that this exponent contributes a suppresion to 
the cross section at leading order, in contrast to the first exponent. This suppresion depends 
crucially on the way the cross section is constructed as discussed in Ref. 0. 

In the third exponent, jg = (as/7r)(3Cir/4) and ■jg = (as/7r)(/3o/4) are anomalous dimen- 
sions of the fields ip for quarks and gluons, respectively. 

The trace in Eq. ( p.7|) is taken in color space. The evolution of the soft function from scale 
Pt/N to pt is given in terms of the matrix Ts, the soft anomalous dimension [^ Q. The symbols 
P and P denote path ordering in the same sense as the variable fi' and against it, respectively. 
The superscript (f) on H, S, and Ts stands for the process fafb -^ JX. The soft anomalous 
dimension matrix, evaluated through the calculation of one-loop eikonal vertex corrections, has 
been presented for all partonic processes in jet production in Ref. [^. It can be written in the 
form 

iXf)KL = iTf,)KL + SKL'^ E %) J [- ln(2z.,) + 1 - TT^] , (2.8) 



TT 
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where the gauge- dependent terms appear exphcitly. We note that all gauge dependence cancels 
out in the cross section. 

In a color basis where the soft anomalous dimension Ts is diagonal, the path-ordered ex- 
ponentials of matrices in Eq. (|2.71) reduce to simple exponentials; however, in practice the 
diagonalization procedure is complicated [p!l|] . 

Note that the Born cross section is given simply by the trace of the lowest-order soft and 
hard functions as 

Ej^^Mi^JX ^ ^rj.^^HUf,^JXgUf,^jx^ (2.9) 

u Pj S 

Since the resummed cross section is given in moment space, one needs to perform an inverse 
transform to calculate the physical cross section. A prescription is then needed to take care of 
the regions of phase space where the running coupling diverges, and the results do depend on 
the prescription. The resummed cross section may also be expanded in powers of a^ and then 
one can perform the inversion order-by-order, thus avoiding the use of a prescription as well 
as the diagonalization procedure mentioned above. We prefer the latter approach, and expand 
the resummed cross section through next-to-next-to-leading order at NLL accuracy. Thus, we 
present the first calculation of NNLO-NLL corrections for jet production. 

In the appendices we present the matrices for the lowest-order hard and soft functions 



(some of which have appeared in Ref. JT^ in the context of dijet rapidity gaps) and the one- 
loop soft anomalous dimensions Ts' f^, ^^, |ll[ for each partonic process in a given color basis. 
The Born cross section, calculated using Eq. ( p. 91) , is also given; it agrees with well-known 



results |T8| , p!9| , pOf , thus providing a check on the lowest-order hard and soft matrices. We then 
calculate the NLO and NNLO threshold corrections to the single-jet inclusive cross section from 
the one- and two-loop expansions of the resummed cross section. In the following section we 
apply those results to high-i^y jet production at the Tevatron. 

3 Jet production at the Tevatron 

In this section the NNLO predictions are compared to data for the single jet inclusive cross 
section at the Tevatron. The NLO predictions have been generated using the EKS jet program 



21| , P^ , P3i with the CTEQ5M parton distributions 0]. The NNLO corrections have been 
calculated in accordance with the expressions contained in Appendices A-E and added to the 
NLO predictions. In all the cases shown here, the renormalization and factorization scales have 
been set equal to each other. 

In Fig. 1 CDF data |jl| are shown compared to the NLO predictions with a scale of Et/2 
in the form of (data-theory) /theory. The excess at high Et relative to the lower Et results is 
clear. The dash-dot curve shows the effect of using 2 E^ as the scale. The scale dependence 
over the E^ range shown varies between about 20% and 28%. Next, the effect of adding in 
the two loop corrections with the scale set to Et/2 is shown by the long dashed line, while the 
corresponding curve for 2Et is shown by the short dashed line. The relative spacing between 
the long and short dashed lines as compared to the solid and dash-dot lines shows that the 
inclusion of the two loop terms has reduced the scale dependence. In addition, the two loop 
corrections are rather small for the scale choice of Et/2. These results are very similar to those 



obtained in Ref. |16] for direct photon production. Using the same notation as in Fig. 1, a 



comparison with D0 jet data [0 is shown in Fig. 2. As is the case with the CDF data in Fig. 
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Figure 1. NLO and NLO + NNLO results for inclusive jet production at the Tevatron compared 
to data from the CDF Collaboration [|l| . 

1, the D0 data show a definite slope with respect to the theory. In fact, allowing for a small 
shift in overall normalization, the deviations of the theory from the data are quite consistent 
for the two data sets. In both cases it is clear that the threshold corrections are not sufficient to 
increase the slope of the theoretical predictions to match that of the data - the high E^ excess 
is still apparent. 

It is normally expected that the effects of the threshold resummation corrections should 
increase with increasing Et since the steeply falling parton distributions restrict the values of 
S4 to be near the parton-parton scattering threshold where the corrections are large. However, 
the results shown in Figs. 1 and 2 show relatively flat corrections. To understand why this 
is so, the ratio of the NNLO correction to the NLO calculation is shown in Fig. 3 for the 
three types of subprocesses - quark-quark (or antiquark), quark-gluon, and gluon-gluon. In 
each case the ratio shows a rise throughout the high E^ region. At high values of E^ the 
corrections to the gluon-gluon subprocess are largest, followed by those for the quark-gluon 
and quark-quark subprocesses. This ordering is in accord with the size of the color factors 
associated with the dominant (In {s/i/p^))/ s^ term. Next, in Fig. 4 the relative contributions 
to the total jet rate from the same three classes of subprocesses are shown. As Et increases, 
the gluon-gluon subprocess falls rapidly, followed by the quark-gluon contributions, so that the 
quark-quark terms dominate at the highest Et- Comparing Figs. 3 and 4, one can see that 
as Et increases, the large gluon-gluon corrections contribute a decreasing share to the total 
NNLO correction while the share from the smaller quark-quark corrections increases, thereby 
yielding the relatively flat behavior shown in Figs. 1 and 2. Predictions for a center-of-mass 
energy of 2 TeV are shown in Fig. 5 where the same behavior is seen to extend past i5T=600 
GeV. We have examined similar predictions for the LHC at 14 TeV and the same pattern is 
seen there, as well. 

In Ref. [Q threshold resummation corrections at leading logarithmic accuracy to dijet 



production in hadronic collisions were considered. There it was found that the the corrections 
to the dijet mass distribution were less than about 5% in the kinematic range considered here 
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Figure 2. NLO and NLO + NNLO results for inclusive jet production at the Tevatron compared 
to data from the D0 Collaboration [H. 
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Figure 3. Ratio of the NNLO corrections to the NLO calculations for various subprocesses. 
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Figure 4. Relative contributions of various subprocesses to the NLO calculation of inclusive jet 
production. 
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Figure 5. Predictions for the NNLO corrections to inclusive jet production at a center-of-mass 
energy of 2 TeV. 



when a scale ^ = Mdijet/'^ was used. This is consistent with our findings, although the details 
of the corrections and some of the techniques used differ in the two cases. 

4 Conclusion 

In this paper the corrections to existing NLO calculations from threshold resummation in jet 
production in hadron-hadron collisions have been studied using the formalism of Refs. [§, ^, 0]. 
Detailed expressions for the NNLO corections to NLL accuracy have been presented for each of 
the parton-parton scattering processes. It is found that adding the NNLO corrections to NLO 
predictions results in a decrease in the scale dependence and that the corrections are rather small 
if the renormalization and factorization scales are chosen to -Et/2. The NNLO corrections are 
smallest for quark-quark subprocesses and increase in size for the quark-gluon and gluon-gluon 
subprocesses. However, the more rapid decrease with increasing Et for the latter subprocesses 
results in a relatively fiat behavior for the NNLO corrections as Et increases. It therefore 
appears that threshold resummation is not sufficient to remove the remaining discrepancies 
between the data and the predictions of QCD for high Et jet production. 
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A NNLO-NLL corrections for quark- ant iquark annihila- 
tion processes 

We begin with the quark- ant iquark annihilation processes, 

q{Pa,ra) + q{Pb^n) ^ q{pi,ri) + q{p2^r2) , (A.l) 

where the pi and r, are momentum and color labels, respectively. We use the notation 

T = In (—\ + TTi , U = \n (—\ + m , (A.2) 



where 



S={Pa+Pbf, t = {Pa-Plf, U={pb-Plf, (A.3) 



are the usual Mandelstam invariants. In the t-channel singlet-octet color basis 

Ci = Clrari^ri,r2 y ^2 = [Tpjr^ra[Tp)r,,r2 y (A. 4) 

where the Tp are the generators of 5*^7(3) in the fundamental representation, the soft matrix 
S at lowest order, with elements 5*^/ = Tr [c^c/], is given by 
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and the one-loop soft anomalous dimension matrix Ts' is |P, P, 25 
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There are three different quark-antiquark subprocesses to consider depending on the quark 
flavors. 



A.l NNLO corrections for QjQj -^ QjQj 

The hard matrix for the process g^gj — *■ qjqj at lowest order, whose elements come from the 



squares of the color-decomposed tree amphtudes, is given by (see also Ref. |17| ) 
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The Born cross section is given, using Eq. ( p.9|) , by 
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The NLO threshold corrections from the one-loop expansion of the resummed cross section 
at NLL accuracy are then 



rf^aH. 



j^ Qili^ljlj 



d?pj 

+ 



^^Jg,^g,5, \ 2Cf 



ln(s4/p 



sa 



1^f\ 3^ (N^ + l), f-t\ (A^f-5). /-u 



-^^-^ns "2^^- 



N. 



^alU_Cl^^(-t_^{t^ + u-^ 



IT Nr 



ln( — 

s 

u\ u^ 



N. 



ln( 

s 



1 

S4 



TT 






Ksa) -Cf 



2^ 3 



Sa 



In ^ + 



in(4Vyi-(f 



.Pt. 



(A.IO) 



From the two-loop expansion of the resummed cross section at NLL accuracy, we obtain the 
following NNLO threshold corrections 
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A. 2 NNLO corrections for QjQj -^ q^qu 
The hard matrix at lowest order is 



ffQjqj^qkqk = q,2 
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The Born cross section is 
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The NNLO corrections are 
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A. 3 NNLO corrections for qjq^ -^ QjQk 
The hard matrix at lowest order is 
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The Born cross section is 
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The NNLO corrections are 
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d^a^^^- - 
"■ ^qjqk^qjqk 

r 

d^pj 

+ 3Cf 



s \ B I 9^2 

^qjqk^qjqk ] ^^F 



ln^(s4/p: 



-2CF\n[^]-\cF- 



S4 

{Nl + 1) 



.Pt. 



Nr 



ln( — 

s 



iN!-5)^^f-u 



N. 



12 



\n\si/pl) 

Sa 
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+ 4C. 



— Cf n — + -Cf + ^^ n — 

2 ls/4 Nr \ s 



(^?-'.„f-uc..n(f; 



iVc V s 



4^-''»'"(S^ 



c. 



2^^111 



VI ^ 



Si 

+ 3Cf + 



-l^MiHS, 



4 
1 

S4 



In^ 



2 / /fF 



■"(i 



1 

Si 



ln(g4/p|) 

S4 



(A.20) 



B NNLO-NLL corrections for quark-quark scattering pro- 
cesses 



Next, we analyze quark-quark scattering processes, 

Q (Pa, To) + q {Pb, n) -^ q (pi, ^i) + q (P2, ^2) • 

In the t-channel singlet-octet color basis 



qq^qq _ (rpc\ (rpc- 

•-l ~ l-^FJnral-f F. 

the soft matrix is given at lowest order by 



r2rb 1 "^2 



aq^qq 



l't' r■^ O'r.i 



rari^rbr2 7 



5 



qq~^qq 



(iVf-l)/4 



(B.l) 
(B.2) 

(B.3) 

(B.4) 
2CfT 

There are two different quark-quark processes to consider depending on the quark flavors. 



N^ 



and the one- loop soft anomalous dimension matrix is P, p5| 



-pqq~*qq _ '^^ 
•s" ~ ~ 



-^{T + U) + 2CfU 2U 



Qftj 



B.l NNLO corrections for QjQj -^ QjQj 



The hard matrix for this process is given at lowest order by 



-"11 -"12 






with 



j^ijij ^ij9j 



'11 



-"12 



Tjqjqj^qjqj 
-"22 



'12 



iVc-- 

tu 



'22 



(B.5) 



2 

iVf 

2Cf 

1^ 



2ciis^ + e) 



1 {s' + e 



NrtU 






m 



u^ 



(B.6) 
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Then, the Born cross section is 



<ij1i-"l3lj 



J^J -T. = O^. 

d-^Pj 
The NLO corrections are 

d^pj 

+ 



Sis,) 



a. 



2ms 



u^ s^ + t^ 



t2 



+ 



U^ 



a 



« ^B 



-a 



n 



ijij^ijij 



ICi 



ln(s4/p 



S4 



-2Cf In 






-2^- 



{Nl 



N^tu 



5(s4) . (B.7) 



N. 



In 



-'VMz^ln 



TT 1 iVc \U 



t\ (s' + e) 



a 



:^^!^,g,^mA'^^ 



-Cf 



s 
.Pt, 



m 

-(^4 



In 



N.. 



+ 
f3o 



+ — In 
2 



-M 






1 

S4 



:b.8) 



The NNLO corrections are 






d^p.j 






In^(g4/PT) 
sa 



+ 3C. 



-2Cf\4^\-\Cf 



Pt 



(iV^ + 1), /-* 



W, 



-w 



12 



ri„2 



In (s^/p. 



S4 



+ 



+ 



a^ 12C| 



ln(- 

u 



2 , f-u\ s 
+ — In - 



■u^ 



«c 



TT 



O" 






^Ct. 



Nr 



-IC..M 



tu 



ln^(s4/p; 



Sa 



+ -Cf + 



{Nl 



N. 



In 



(^.(^)...J^^ 



.Pt. 



\n[^]+lcFpM^ 



.Pt. 



.Pt. 



+ 



ai 16C| 



TT 



+ 



2 iVc 



ln(- 



t\ is^ + e] 



w 



Nr ^\ S ) tu 



"■(f 



ln(s4/p^ 

S4 



ln(s4/p^ 
Sa 



TT 



O" 



.B 



Cf 



2CF\n(Pf\+3CF + ^ 



In 



2 / A^F 



-|^^^»'"(|)Kf 



1 

S4 



:b.9) 



B.2 NNLO corrections for qjQk -^ Qjqk 
The hard matrix at lowest order is 

~ 2{s^ + u^)/{N^t^) 



JJQjQk~*qiQk 



a. 











(B.IO) 
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The Born cross section is then given by 



The NLO corrections are 
d?p^ 



'qjqk 



-af^. 



QjQk- 



^q^qk^i'^) = ^s /^2/ t2 ^ (^4) 



+ 



qjqk^Qjqk 1 ^ 



TT 9j9fc 



ln(s4/pg 
S4 



-2C..„f§ 



-«(S4) 



The NNLO corrections 

-'qjqk 









04 j + 






pfj-f'"(S, 



^ '^ %gfc- 



rf^p^ 



»j 



are 



O", 



"JiSfc- 



-<?i9fc ] 2C^ 



In^ 



(g4/PT) 

S4 



-3Ci7 



-^•"(t)-^-^-T'"(T; 



+ 4Ci. 



c-3)^ /-^\ /^o' 
A/'e V s / 12 

/^2 \ Q / Ar2 



, {Nl - 3) 



In^ 



(g4/PT^ 



S4 






ln(s4/p^ 

S4 



V^Ty 



■"(s: 



4) 



+ 






ln^f4 



PTy 



1' 

S4. 



^)ln(^ 



1 

.S4 



fB.lll 



1 ■ 

S4- 



ln(s4/p^ 

Si 



(B.12) 



(B.13) 



C NNLO-NLL corrections for qq -^ gg and gg -^ qq 



Next, we analyze the processes qq -^ gg and gg -^ qq. 
For the process 

q (Pa, Ta) + q (Pb, Tb) ^ g (pi, Ti) + g (p2, ^2) 



in the s-channel color basis 



qq^gg _ x x < 

•^1 "rari,"rir2 ? "^l 



.qq^gg 



rf"i"2"(T;; 



rtra ' 



^qq^gg _ j^sr^r^c/rpc-. 



nra ' 



(C.l) 
(C.2) 
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where d"*'^ and f"''"^ are the totally symmetric and antisymmetric SU{3) invariant tensors, the 
soft matrix at lowest order is 



(^QQ^gg 



Nc{Nl - 1) 

{Nl-A){Nl-l)/{2N,) 









Nc{Nl - l)/2 



and the one-loop soft anomalous dimension matrix is 





J- s' 



TT 










U -T 



2{U-T) ^{U-T) ^{T+U) 
These same two matrices also describe the time-reversed process P, ^ 

giPuri) + gip2,r2) -^ qiPa,ra) + q{,Ph,n) . 



(C.3) 



(C.4) 



(C.5) 



C.l NNLO corrections for qq -^ gg 

For this process the hard matrix at lowest order is given by 



H^^gg = ^l 



mi 



qq^gg 



TTqq--*gg 

-"12 

Tjqq^gg 

-"13 



Tjqq^gg 

-"12 

Tjqq^gg 

-"22 

Tjqq^gg 

-"23 



zjqq^gg 

-"13 

rrqq^gg 

-"23 

TTqq^gg 
-"33 



with 



j^qq^gg 



ni 



2M 



u t 
t u 



Tjqq^gg 



'12 

rgq- 

'22 



Tjqq^gg 






Tjqq^gg 

-"13 



rrqq^gg 

-"23 

zrqq^gg 

-"33 



iu^ - 1^) 



[u 



t) 



2iV3 



tu 



AT rrqq^gg 
^^c-"l3 1 



2NI tu 



s^ 4 tu 
+ 



ms^ 



ATS 



ATS 

c 



The Born cross section is given by 

,2(A^c'-l) 



E 



d^a^ 



j- 



qq^gg 



a. 



B_ 

qq^gg 



E, 



(Pp.] 
The NLO corrections are 

'-*' ^ an — ¥nn "-*• 



5{s,) 



a. 



NrS 



{Ni - 1) {e + u^ 



2m 



tu 



qq^gg 



* _B 



c^^Pj 



-o". 



TT 



qq^gg 



{ACf - 2Ca) 



ln(s4/p|) 



S4 



S2 



(C.6) 



(C.7) 



5{s^) . (C.8) 
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+ ^/ (^ 



2 



^ Os^B 






The NNLO corrections are 



tu \sj 2 \ 



- 1) ({v?-f) ^ 20 



1 

.S4- 

tu 
2 Ip^, 



"^)'"(i; 






(i^p„ 



'J 



- (?)^.- 






*99 



^(4Cp-: 



'In^l 



Sa 



+ 



3(2Ci 



■> — Ca) 






[sj/pV ^ 



>t) 



+ 



+ A 



'0 



9 N 



rin^i 



S4 



« 



+ ^ 3 (2Cp - Ca 

Vf-l 



{Nl 



{Ni - 1) {e + -2 

2iV2 



tM V s ; 



iu' ^2 



tu 



sa 



In^ 



as\2 



+ (- 






fc. - ^^ 



.M'"^^- 



a 



H — - 



.^^'^^'"(S! 



;|),34c.-^),„(|)}[!^ 

2iV2 



lnf4^ 



4) 



+ 



(Al 



±^lnp 

t- 



'l)( {u'-t^) 20 



t-u 



"^)in(^; 



ln(s4/p^ 

S4 



+ 



;^)VS._ /c.U.„(ii).3c. + f]>.^(f 



^Crf„.ng).nf4 



vPt/ 



1 
S4J 



1 ■ 

S4- 



(C.9) 



(CIO) 



C.2 NNLO corrections for gg -^ qq 

Note that apart from color factors, the hard matrix for this process is the same as for qq -^ gg: 



H99- 



■ qq 



(iV2 - 1)2 



H 



qq^99 



(C.ll) 



Then for the process 



g {Pa, r-a) + g (Pb, n) -> g (pun) + q (^2, ^2) 



(C.12) 
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the Born cross section (symmetric in t,u) is given 






by 



+££ 



— ^B 



■- a 



The NLO corrections 

99^11 s B 



{Nl 



P,-l)s 



\Nl-l)it' + 



u' 



2iV2 tu 



ill' 
1 + 2- 5{s^). (C.13) 



B 

^ 99^11 



are 
^(4Ca 



-2Cf 



ln(s4/p^ 



4) 



S4 



+ 



|c.-(2C.-C.)ln(li)-2C.l„(|)" 






a|r_(iy^-2)(t^+^j /p|\ (Ai_ 

TT \ 2Ar2 ^y (^ g j 2 



■ 1 

.54- 



- 1) ({u^-e) ^ 20 



+ 

2(u 



tu 



T^X, 



'^<^J(s.)\-C^\4'^ 



The NNLO corrections are 



sA$)Mi. 



-1 "• ^99^<1<1 



d?pj 



= f^V 



TT / ^^ 



^A\{ACA-2CFf 



In^l 



>4/Pq 



'tJ 



+ 



2>(2Ca - C, 



3 



S4 

•.-c.).,>(?i)-2c..,>(|; 



+ /9o(-CA + ^Ci. 



In^ 



(S4/P?^1 



f^TJ 



S4 



7r2 



5 (2Ca - 



^i - 1; 



2iV2 






(Ag 



(w 



s 



In^ 



{sa/pI 



>t) 



S4 



+ ^ 



7; ^fs-??]^^^^ 



\ /I / I \^ / \ o / ^ 

- P,{2Ca - Cf) In (41 + ?/5o(2C^ - C^) In (4 

\prp J Z \PT y 



ln(s4/p|, 

S4 



1- 

.S4-- 



(C.14) 






0' 



|4C^lnf4 






(u^- ^2 



^Vnf-^ 



2 y t-u s J \t , 



ln(g4/pj 
S4 



4) 



1 
.S4J 



+ 



(C.15) 
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D NNLO-NLL corrections for qg -^ qg 

Here we discuss quark-gluon scattering, 

Q (Pa, ra) + g (Pb, n) -^ q (pi, ri) + g (p2, ^2) • 
In the t-channel color basis 



qg^gg _ x x 

"--1 "rari"rf,r2 5 

the soft matrix at lowest order is 

N,{N^ - 1 



(ig^ig 
^2 



dnr,c^j.c^^^ 



ra ' 



O 



.Qg^qg 



^r'-nnx.ra 



gqg^qg 








{N^-A){N!-l)/{2N^ 













Nc{N! - l)/2 



and the one-loop soft anomalous dimension matrix is 

{Cf + Ca)T 



rqg^qg _ '^^ 

s' — ~ 

71 





2U 



CpT + ^U 






U 

Cau 
2 ^ 



2Nc 



CpT+^^U 



The hard matrix at lowest order is given by (see also Ref. ||T7[| ) 



jjqg-^qg = ^2 



Tjqg^qg rjqg^qg rjqg^qg 



'12 



'13 



jjqg^qg j^qg^qg j^qg^qg 



'12 



'22 



'23 



with 



rrqg^qg 
-"11 



Tjqg^qg 



Tjqg^qg rjqg^qg rrqg^qg 
-"13 -"23 -"33 



1 



'12 

'22 



TTqg~*qg 



2A/'3(A/'2 - 1) \su 



Tjqg^qg 

-"13 



Tjqg^qg 

-"23 

zrqg^qg 

-^33 



Nl{N2 - 1) 
qg- 

13 

1 



2s u s ' 
^T^2s~2u. 



AT TJ-qg^qg 
^^c-"l3 ) 



Nc{N^ - 1) 



Asu 



f 



t2 2su 



The Born cross section is 



E 



dS^B 



J- 



qg^qg 



d?p 



a 



J 



B 

qg~^qg 



5{s^) 



2I 

OiZ- 



^ m 



2m SU t2 



Ksa) ■ 



(D.l) 



(D.2) 



(D.3) 



[D.4) 



(D.5) 



[D.6) 



(D.7) 
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The NLO corrections 

(Ppj 



are 






'":^<^,AiCF + c.) 



In(s4/pj 

Si 



4) 



+ 



+ 



-{CF + CA)ln(^ 






+ — 

TT 



1 









vPt 



' 4 



t- 



t 



+ 



u 



Ca 



xJ^' 



kn^ \ 






-2 



+ 



a. 



The NNLO corrections are 

gg- 
d^pj 



c.i„(^; 



-C.lu(^)+^C.+ 



Is 2uJ \ s 
lnf4 



1 

-S4- 



1 1 

LS4. 



:f)4-(f: 



(D.8) 






= r^v 



, TT 



0-, 



"JS- 



.qg{liCF + CA? 



In^l 



!g4/PT' 



S4 



+ -{Cf + Ca) 



-(C. + C.)ln(|)+C.(-21n(^; 



4y 



-2C.lnf^)-^l[i^^!i^ 
V s / 3 S4 



TT"' 2 



P^) 



1 /t^ 

c 



a 



■A 



In r^^ 



S4 



;w-) 



2s ~ 2m 
'In^i 



Hr, 



'mIpI 



>t) 



S4 






. TT 



<^JCf + Ca) I (4C^ - Ca) In (^) + {ACa - C^) In (-^) +\c^ + ^ 






'<t)n^»'"(s: 



ln(g4/pj 



4) 



at 

H — ^ 
2 



-2)(C,. + Ca)' ^ 



+ (c.in(^; 



Ca 






2s u 
-1- — + 



t 2s 2m 



llnf^^ 



InC^ 

S 



vW-^J 



su 



-2 +2 1 



^?)1WI)[^ 












■i-KDKi; 



■ 1 
.S4J 



+ 



(D.9) 
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E NNLO-NLL corrections for gg -^ gg 

Finally, we consider gluon-gluon scattering, 

g {pa, Ta) + g (Pb, n) -^ g (pi, n) + g (p2, ^2) . 



(E.l) 



The color decomposition for this process is by far the most complicated. For simplicity in this 
section we use A*",, = 3 explicitly. A complete color basis for the process gg -^ gg is given by 
the eight color structures [§ 









fraTbl jrir2l _ ffaTbl £rir2l 



^2 



c: 



99^99 



-eraTtl jrir2l i ^farbl J!rir2l 



0rar\0rbr2 i 



39^99 



C 



„99^99 



J99^99 



_jraricrrbr2C 99^99 . 

5 

1 



_ fraric erbT2C 
O 

_ _ fTaVlC £ri,r2C 

O 



l"rari,0rir2 + "rar2'^ri,ri j o '""''i '^6^2 



_jraricrri,r2C 



(E.2) 



The soft matrix at lowest order is 

g99^99 



099^99 n 
•^3x3 U3x5 



J5x3 



099-^99 



'5x5 



where 



q99~*99 



■^3x3 



5 





■ 





5 











5 _ 



q99~*99 



1 '-'5x5 



10 

8 

8 
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The one-loop soft anomalous dimension matrix is 



-p99^9g 



r99^99 n 

■■-3x3 U 



-"5x3 



3x5 
-p99^99 

'- 5x5 



with 



and 



-p99-*99 
'- 3x3 



TT 



3T 

3U 

3{T + U) 



r99-'99 

'- 5x5 



a.9 



vr 



6T 

3T + 3^ 



-6f/ 



3U_ 

4 








_ —EL 

2 2 

■f 3T + f 

■^ 

-^ 

^ 3 





-3f/ 



3f/ 

4(7 








9C/ 

4 

5 



-2T + AU 



(E.3) 



(E.4) 



(E.5) 



(E.6) 



(E.7) 



The hard matrix at lowest order is |17 



where 



with 



a„ 



^99^99 



'5x5 



tt99^99 

-nil 

tt99^99 
-"12 

TTg9~*99 
-"13 



tt99^99 
-"15 



-"3x3 



-"5x3 



H, 



-"3x5 
99^99 



Jj99^99 



'12 



'99- 
22 


*99 


23 


*99 







•99- 
25 


*99 



5x5 

tt99^99 
-"13 

tt99~^99 
-"23 

tt99-^99 
-"33 



tt99^99 
-"35 



i^ffr^^ 



15 



i/f^^^^ 



'25 



/7|3^^^ 



'35 



/7f ^^^ 



'55 



22 



(E.8) 



(E.9) 



i7 



99-*99 
11 



tt99-^99 
-"12 



tt99^99 
-"13 
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tw st 



16 y s^ M^ swy 
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2^11 ' 

9 /st tM M^ S^^ 

^ Ht-^ + — - — 



32 \m2 



st tu 
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if. 
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22 
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tt99^99 
-"25 



6 

st 



F^ff" 



►99 



27 9 /su tw 

64 ~ 16 V^ ^ 4? ^ 4^2 



_9_/n2 
^ 32 \ st 



s 
tw 



2SM , 



tt99^99 
-"35 



TT99^g9 
"3^13 , 



TT9g-*99 
-"55 



9 ^^ 



The Born cross section is 
Ej 



(Pa^ 



99^99 — ^B 

~ 99-*99 



E, 



The NLO corrections are 



Ksa) 



QJo 



27 9 f su tu st\ 
+ - + - 
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ln(s4/p^) 
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TT I 8 
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-t 
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2 V s 



27 



su 



u^ 



tu 



su 



4 V s 



5(S4) 



(E.IO) 



(E.ll) 
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(E.12) 
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The NNLO corrections are 






aa^gg 



(Pp^ 



'J 



-^'. 



0-, 



'gg^aa | ^^A 



"In^l 



>4/Pt! 



S4 
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^M'l 
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+ ' ' 



2 
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.B 
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-99 \ ^A 
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t-u st 
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su 
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S4 
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2 Ui ^ 



-M^ 
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■u- 
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27 
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In I 



-M^ 



+ 






^^9-99 W^ 



,v)][t-H?- 



27 /SM tn st \ 9 /u 



2 1st ^ 



e 



tu 2su 



'1] In f /^' 
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st 

1- 

■S4J + 



S4 

In(s4/Pr) 
t-u u"^ i 
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st tu 



ln(s4/p^ 
S4 



4)1 



(E.13) 
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